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We present approximate methods for calculating the three-dimensional electric potentials of finite 
surface electrodes including gaps between electrodes, and estimate the effects of finite electrode 
thickness and an underlying dielectric substrate. As an example we optimize a radio-frequency 
surface-electrode ring ion trap, and find that each of these factors reduces the trapping secular fre- 
quencies by less than 5% in realistic situations. This small magnitude validates the usual assumption 
of neglecting the influences of gaps between electrodes and finite electrode extent. 
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Trapped ions have shown to be excellent candidates 
for building quantum information processing experi- 
ments [IH!] • It is possible to trap individual ions [6] , 
manipulate their motional and internal degrees of free- 
dom with outstanding fidelity [B], and read out their 
full quantum- mechanical state [7]. In order to construct 
small and scalable setups of trapped ions it is necessary to 
build and improve radio-frequency (rf) Paul microtraps 
constructed from surface electrodes [SJ |U] • The design of 
surface electrodes is well understood [TUWT^] as long as 
several simplifying assumptions are made. In this work 
we discuss two of these assumptions, present techniques 
for more accurate calculations leading to better designs, 
and give estimates of how good these assumptions are in 
practice. 



The particular assumptions under study in this work 
are (i) that there are no finite gaps between surface elec- 
trodes and (ii) that the surface electrodes extend to in- 
finity and cover the entire plane. Such planar infinite 
gapless electrodes are assumed in order to simplify theo- 
retical work [TUHT^] and to achieve analytic and general 
results. By contrast, experimental setups are often simu- 
lated with boundary-clement (BEM) or similar numerical 
methods |13[ 1 14] , which are much more demanding and 
only give narrow results for specific setups. All of the re- 
sults achieved in this work can be further improved upon 
by such numerical methods; our goal here is merely to 
estimate the size of the corrections due to the aforemen- 
tioned assumptions and to suggest fast approximations 
for taking the leading corrections into account during the 
design phase of an ion trap setup. 



This article is organized as follows. In Section |T] we 
summarize the tools used for calculating the electrostat- 
ics of planar electrodes. Section [TT] looks at the influence 
of gaps between electrodes, and Section |III| at the influ- 
ence of the finiteness of the surrounding electrode. 
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I. SURFACE CHARGE DENSITY 

The mathematical tools for calculating the electric 
fields of planar infinite gapless electrodes have been pre- 
sented in detail elsewhere [TUl US]- They assume that 
there is an infinite plane in space, here the xy plane, 
on which the electric potential $(r) is fully known since 
it is determined from the shapes of the gapless surface 
electrodes which are individually connected to externally 
generated voltage sources (dc or rf). This potential is 
then mathematically extended away from the electrode 
plane through Laplace's equation V 2 <I>(r) = satisfied in 
vacuum. In all generality this extension can be expressed 
as a surface integral 

$(x,y,z) = <f>(x',y')G(x-x',y-y',z)dx'dy', (1) 
Jr 2 

with the surface potential <f>(x,y) = &(x,y,0) and the 
Green's function 



G(x,y,z) 



2ir(x 2 + y 2 + z 2 ) 3 / 2 ' 



(2) 



Equation ([I]) is sufficient for performing calculations un- 
der the approximations of infinite electrodes and no gaps, 
since the integral is completely determined through the 
full knowledge of 4>{x, y) and thus the potential $(f) can 
be calculated anywhere in space. Many other tools devel- 
oped for such planar infinite gapless electrodes are deriva- 
tions of Eq. Q. 

Both approximations studied here, gaps between elec- 
trodes and empty space surrounding a finite electrode, 
have in common that their relaxation leads to areas in the 
xy plane where the electric potential <fr(x, y) is not known 
a priori and hence Eq. ([!]) can no longer be used with- 
out first determining the surface potential in the open 
areas. Here we derive additional constraints on the po- 
tential $(f) which are still satisfied even if such open 
(electrode-free) areas are present in the xy plane. 

If we view the planar electrodes as an infinitely thin 
sheet floating in vacuum at z = 0, Eq. ([!]) yields a sym- 
metric potential $>(x,y, — z) = <f>(x,y, z). We define the 
surface charge density of the sheet [TB] 



<j(x,y) 



lim 

z ->o+ 



d$(x,y,z) 
dz 



(3) 
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where £o is the permittivity of free space. As long as the 
assumptions of no gaps and infinite electrodes are satis- 
fied, a(x, y) can take arbitrary values on the electrodes 
since charge carriers in metallic electrodes will move to 
wherever necessary for generating the surface potential 
4>{x, y). However, whenever there is a hole in the surface 
electrodes, for instance within a gap between two surface 
electrodes, a(x,y) must be zero in this hole since charge 
carriers cannot exist in the vacuum outside of metallic 
electrodes. This condition gives us a way of calculating 
the three-dimensional electric potential of gapped and fi- 
nite planar electrodes. It is important to note that the 
surface charge density ([3| and the condition that it van- 
ish outside of the physical electrodes are only valid when 
the electrodes are infinitely thin and the potential sym- 
metry $(:£, y, —z) = <5>(x, y, z) is satisfied. In Section HE 



we study the effect of thick electrodes as well as substrate 
dielectrics on the conclusions reached from thin-electrode 
calculations. 

The surface charge density (constrained to zero out- 
side of the surface electrodes) and the surface electric 
potential (constrained to given values inside of the sur- 
face electrodes) are related through 



(x,y) = 4> + 



1 



4tt£ 



<j(x', y')dx'dy' 



£o 



cr(x, y) = — - hm 

7T z-S-0 



yj{x-x') 2 + {y-y') 2 
{x',y')&x'dy' 



(x - x') 2 + (y - y') 2 - 2z 2 
[(x-x') 2 + {y-y') 2 + z 2 } 5 / 2 ' 



(4a) 



(4b) 



The order of the limit and integral in Eq. ( 4b ) can in gen- 



eral not be interchanged. In order to compute the surface 
charge density from the surface potential, it is necessary 
to compute a related quantity (the z component of the 
electric field) for z ^ and then take the limit z ^ 0. 
This makes many calculations involving gapped or finite 
electrodes much more demanding than they first appear. 
In practice an approximation is often sufficient: defin- 
ing a small region A(x,y) in the xy plane close to the 
point (x,y), we can split the integral into two parts, in- 
tegrating over A(x, y) and M 2 \ A(x, y) separately. In the 
latter part the limit and integral may be interchanged, 
thus eliminating the need for a three-dimensional inter- 
mediate; in the former part a series expansion of 4>{x' , y') 
around the point {x,y) is usually sufficient. If A(x,y) is 
a square of sides 2d x 2d centered at (x, y), then the local 
(singular) contribution to Eq. (4b) is found from 



lim 



rx-\-d rV+d 

\ Ax' / &y'<t>{x\y' 

Jx—d Jy—d 



Ay/2 



(x - x 1 ) 2 + (y - y') 2 - 2z 2 
\{x - x') 2 + (y - y') 2 + z 2 p/2 

{x,y) + 2dsmh- 1 (l)\7 2 c/)(x, y) + 0{d 3 ). (5) 



As d — ¥ the first term of the series in Eq. |5]) diverges 




FIG. 1: Perspective view of an infinitely long straight gap of 
width g between two semi-infinite planar electrodes. The left 
(gray) electrode at x < —g/2 is grounded; the right (green) 
electrode at x > g/2 is held at unit potential. 



but is cancelled out by an equal and opposite divergent 
term in the integral over M 2 \ A(x,y). 



II. GAPS BETWEEN ELECTRODES 

Analytic studies of ion traps built with surface elec- 
trodes are often simplified to gapless electrodes [T0HI2] . 
As a first improvement over this approximation Ref. 
introduced an interpolation method for calculating the 
influence of gaps between electrodes to lowest order. We 
show here that this method fails to take into account 
the dominant "gap polarization" induced by surround- 
ing electrodes, and propose an improved method. 



A. infinite straight gap 

We begin our discussion of gaps between surface elec- 
trodes with the simplest possible case: an infinite straight 
gap between two thin infinite semi-planes at different po- 
tentials (Fig.[T]). Its scaled surface potential 



if x < -g/2 



if»(x) = 



isin" 1 ^ ii\x\<g/2 (6) 



9/2 



if x > g/2 



is plotted in Fig. [2] Its associated surface charge density 
is calculated from Eq. (4b I, 



4e 



- gap (x) - 



) 



if x < -g/2 

if |x| < g/2 
if x > g/2. 



(7) 



FIG. 2: Scaled surface electric potentials in the setup of Fig. [T] 
Solid line: gap interpolation potential for thin electrodes, 
Eq. Q. Dashed line: gap polarization potential, Eq. Q. 



The surface potential ^ matches the desired electrode 
potentials on the electrodes (|x| > g/2), and the surface 
charge density |7| vanishes in the gap (|x| < g/2) as re- 
quired. In fact Eq. ([6| gives the only potential which sat- 
isfies both of these conditions, and it coincides with the 
gap potential as determined by numerical finite-element 
calculations. However, in what follows we will see that 
this potential is not sufficient to describe arbitrary gaps. 

As Eq. ^ gives the surface potential in the entire xy 
plane, we can now use Eq. ([I]) to calculate its associated 
three-dimensional electric potential 



*r p (o = 2 



i i 




n 



g 2 sin(2i?) 
16nr 2 



+ 0[(g/r)% (8) 



with 



■ sin i?, 



= rcos$, and Z 



_ \z\+ix _ re^_ 
9/2 3/2 ■ 

The first two terms in the series expansion are due to the 
infinite electrodes, and the remainder to the finite gap. 
In the same geometry the surface electric potential 




leads to the surface charge density 



< ol (z) 




a 1 



2N 



if \x\ < g/2 
if \x\ > g/2 

if \x\ < g/2 
if |sc| > g/2. 



(9) 



(10) 



While the isolated straight gap of Fig. [T] is fully described 
by Eq. Q, the "polarization potential" of Eq. Q will 
appear in more complicated electrode setups. It vanishes 
on the electrodes (\x\ > g/2), and the associated charge 
density is constant inside the gap (|x| < g/2). Whenever 
an unphysical charge density is induced within a gap, 



FIG. 3: Perspective view of an infinitely long strip electrode of 
width w at unit potential (green) within an infinite grounded 
plane (gray), separated by gaps of width g. The strip width 
w is measured between gap centers. 



usually by nearby electrodes via Eq. (4b), we can add 
a surface potential inspired by Eq. |9]) in order to bring 
the in-gap charge density back to zero without modifying 
the potential on the electrodes. The associated three- 
dimensional electric potential is 



$P ol (f) =5R(Vl + Z 2 -Z 



(/COS1? 

4r 



0[(g/r) 



where 5ft refers to the real part. Comparing the far-field 
expansions in Eq. ^ and Eq. (Ill we already see that 
whenever a gap polarization potential similar to Eq. |9]) 
is present its far field will dominate over that of the gap 
interpolation potential (|6| because of its slower decay 
with r. 

The two gap potentials Q and ^ will be shown to 
suffice for describing thin gaps by comparison to numer- 
ical calculations (see Fig. [7]) . We will therefore not need 
any further potentials in the following sections. 



B. two parallel straight gaps 

As an application of the surface potentials (|6| and ^ 
we study a strip electrode of width w at unit potential 
separated from an infinite grounded plane by two gaps of 
width g <C w (Fig. [3]) . The surface potential is approxi- 
mately 



tf»(x 



+ a x 



+ o)-0f P (£-77) 



-7,) + <f> P q °\x--) , (12) 



where a is a parameter to be determined. It is not a 
priori clear that a = as was the case for the isolated 
gap in Section II A for electrodes differing from the setup 



4 



of Fig. [T] the potential at the center of the gap need not 
necessarily be the average of the two adjoining electrode 
potentials. Electrodes which are farther away, as well 
as the curvature of the gap (see Section II C), can in- 
duce such a polarization potential since there are no free 
charges within the gap to counteract these influences. We 
determine the correct value for a by setting the surface 
charge densities at the centers of the two gaps to zero, as 
required physically: 



strip 

w,9 



( X = ± — 1 



4£o 

- 



1 — 2i:aw/g 
\J 4w 2 — g 2 



2-KOt 

9 



0, (13) 



giving 



2n(w - y/Aw 2 - g 2 ) 



9 
2nw 



0[{g/wf]. (14) 



In Fig. [7] we will show that this potential matches that 
determined by numerical simulations for g <C w. 

The gap contribution to the electric potential far from 
the strip electrode is dominated by the gap polarization 
component, 



<i p M)= 1 



irag\ wcos$ 
~2w) 



nr 
1 



+ 0[(w/r) 3 } 

g 2 \ w cos $ 
Aw 2 



(15) 



This setup demonstrates that when taking gaps between 
electrodes into account, it is crucial to calculate the gap 



polarization since its contribution to the electric poten- 
tial is generally more important than that of the gap 
interpolation potential ([6]). Furthermore the gap- width 
dependent prefactor of 1 — g 2 /(Aw 2 ) in Eq. (15) gives a 
first quantitative idea of how important gap potentials 
are. Even for g = w/2, where the width of the gaps is 
equal to the width of the strip electrode, its value is still 
94%. 

We have limited this discussion to nulling the charge 
density at the center of the gap. For increased accu- 
racy one can derive generalizations of the polarization 
potential with in-gap charge distributions proportional 
to successively higher powers of x which can serve to null 
successively higher moments of the in-gap charge. How- 
ever, the electric potentials associated with such higher- 
order gap potentials drop off faster with distance from 
the gap than Eq. ( 11 ) and will therefore be neglected in 
this work. 



C. ring ion trap with gaps 



As a more practical application of gap potentials we 
study an optimized radio- frequency ring ion trap [TUl [H] 
in the presence of gaps between ground and rf electrodes 
(Fig. [4]) . Assuming thin gaps and inspired by the results 
of Section [lIB| we take the cylindrically symmetric scaled 
in-plane potential as 



'0 

6 ring (r) = I 1 




,ga P(r 



e p (i? 2 



if r < R ± - § 
i? 1 ) + a 1 0P° 1 (r-i? 1 ) if|r-fli|< 



9 

APol 



if #1 + f < r < R 2 - f 



(16) 



r) + a 2 0P ol (r - R 2 ) ii\r-R 2 \<\ 
if r > R 2 + f. 



The surface charge densities at the gap centers are, to 
leading orders in the gap width g, 



a™z{Ri) w -2e 



K 



( 4fli_R 2 
\{Ri+R2)' 



^ UR1+R2) 2 



tt(R 1 +R 2 ) ^{R 2 -R l ) 



4R 1 R 2 
(R2-R1)' 



2^R~i g~ + 2(R 1 +R 2 ) 2 {R 2 -R 1 



(17a) 



a lins (R 2 ) w -2s 



K 1 ARtR2 
' V 1 (R1+R2) 



p 1 4R t R 2 
U \ (R1 + R2) 



+ 



n(Rx + R 2 ) n(R 2 -Ri) 

l0 §(32k) 2a 2 , 9RiaiE(-^^ 
2irR 2 g 



2(R 1 +R 2 ) 2 {R 2 -R 1 ) 
(17b) 



where K(m) and E(m) are complete elliptic integrals 
of the first and second kind, respectively. Setting 



KRi 



ring 



(R 2 



allows us to compute the 



amplitudes ct\ and a 2 necessary for nulling the surface 
charge densities at the centers of the two gaps, to lead- 
ing orders in g. 



We do not have a closed expression for the 3D potential 
due to Eq. (16); but on the z axis (r = 0) it takes the 
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FIG. 4: Perspective view of a planar ring ion (Paul) trap 
with gaps. An rf electrode (green) is embedded in an infinite 
grounded plane (gray) and surrounded by gaps of width g. 
The radii Ri t 2 are measured to the gap centers. 
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FIG. 5: Optimized planar ring trap with gapped electrodes, 
trapping at height z above the electrode plane. As a function 
of gap width g, the optimal radii and the maximal achievable 
curvature k differ from the idealized result at g = 0. 



simple form (to leading orders in g) 



$""g(r = 0,2) = 



.9 



4 

2 'z| 



16 



(i? 2 +Z 2 ) 3 / 2 

2R\ - z 2 



^/R 2 2 + z 

(R 2 + Z 2 ) 3 / 2 



2R 



^2 1 



{R\ + Z 2 fl 2 (i?2 + z 2)5/2 



(18) 



Similar to the system in Section [II B| for small gap widths 
g the influence of the gap polarization potentials [second 
term in Eq. ( IS )] is found to dominate over that of the 
gap interpolation potential [third term in Eq. (18)]. 

Eq. ( 18 ) allows us to calculate the optimal rf ring radii 
as a function of gap width. The ponderomotive pseu- 



dopotential of an ion of mass m and charge q moving in 
an electric field $(r,i) = U l{ cos(fi rf i) x $scaled^ ig jg] 



*(r) 



scaled 



4m£! 2 f 



(19) 



For a given gap width g and desired trapping height z wc 
determine the radii R 12 for which \&(r = Q,z) = and 
which maximize the dimensionless trap curvature |12) 



k = z 2 |detif$ 



scaled 



(r) 



1/3 
trap 



(20) 



defined through the Hessian determinant of the scaled 
potential at the trapping site. Fig. [5] shows these quan- 
tities relative to their well-known values at the gapless 
limit (? = 0. We find that for finite gap widths g the in- 
ner radius R\ must be slightly increased; but the limiting 
value R\{g = 0) = 0.678z always lies within the gap. The 
same is true for the outer radius, whose optimized value 
is even less sensitive to the presence of a gap. Finally, the 
trap curvature k is bound to decrease from its gap-free 
value, but this decrease is limited to a few percent for 
reasonable gap widths. 



D. arbitrary gap shapes 

For arbitrary gap shapes, as they naturally appear e.g. 
in surface electrode optimizations |12j . simple homoge- 



neous parametrizations such as Eqs. (12) and (16) are 
not applicable. Instead, as long as the gaps are suffi- 
ciently smooth and have radii of curvature much larger 
than the gap width, we propose the following scheme. 
The center of each gap is parametrized by a curve 7*(i) 
winding counterclockwise around an rf electrode, and a 
normalized direction e(t) — f"^ 1 ] • 7 / (*)/||7 / (i)|| is de- 
fined at each point perpendicular to the gap direction. 
Each point p within the gap can now be described by lo- 
cal coordinates: p(t,u) = j(t)+ue(t). Inspired by the re- 
sults of the preceding sections we propose to parametrize 
the scaled gap potential as 

^(t, u) = 4%*(u) + a(t) x 0Pgj( u ) (21) 

whenever |u| < g(t)/2. As long as a(t) is a sufficiently 
smooth function of t, with variations on length scales 
much larger than the gap width g(t), this potential will 
be a good representation of the true gap potential as 
determined by numerical inversion techniques. The full 
in-plane potential is given by a term like Eq. ( 21 ) for each 



gap curve plus the potential on the electrodes. The exact 
way in which gap intersections and forks are incorporated 
into the above parametrization is largely irrelevant, since 
they represent zero-dimensional objects in the electrode 
plane and their 3D electrostatic potentials thus decrease 
faster with distance than those of the corresponding one- 
dimensional gap potentials, Eqs. ^ and (11). 



As can be seen from comparing Eqs. ( 13 ) and ( 17 ), the 



surface charge density at the gap center depends strongly 
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dielectric substrate 



FIG. 6: Cross section through a gap of width g between two 
coplanar electrodes of thickness h mounted on a substrate of 
dielectric constant e s . The gap susceptibilities of Fig. [7] are 
computed in the top plane indicated by a dotted line. 



on the gap shape and curvature. For general setups the 
surface charge density at the gap centers must be calcu- 
lated using an approximation such as Eq. (5]), giving a 
functional dependence of the a centel (t) on the functions 
a(t). A numerical inversion then yields the necessary 
a(t) for nulling the surface charge densities at all gap 
centers. Since this inversion procedure involves only one- 
dimensional gaps it is expected to be much faster than a 
full BEM calculation. 



E. electrode thickness and substrate permittivity 

One of the main assumptions of this section is that 
the surface electrodes are infinitely thin and lie in one 
plane. This assumption allowed us to define a surface 
charge density <j(x, y) in the electrode plane, which could 
then be set to zero outside of the metallic electrodes in 
order to determine the surface potentials. In realistic se- 
tups, however, electrodes have a finite thickness and are 
often mounted on a dielectric substrate. Under these cir- 
cumstances numerical inversion techniques are necessary 
to determine exact electric potentials. But an approx- 
imate method can nonetheless be derived in the spirit 
of the preceding sections. In Eqs. (151 and (18 1 we had 



found that the far-field of the electrode gaps is domi- 
nated by the "gap polarization" , whose amplitude a must 
in general be determined by an inversion procedure (Sec- 
tion II D ) . Here we suggest that to leading order the cross 
sectional shape of the gap (Fig. |6| merely introduces an 
effective gap susceptibility which scales the gap polariza- 
tion potential by a numerical prefactor. 

We have used two finite-element computer codes 
(translation invariant and cylindrically symmetric) to 
calculate the effects of finite gap aspect ratios and sub- 
strate dielectric constants (Fig. [6| on the surface elec- 
tric potential (along the dotted line in Fig. [6j. As in 
thin-electrode calculations, this surface potential in the 
gap can be decomposed into a gap interpolation com- 




0.01 0.1 1 

gap aspect ratio h/g 



FIG. 7: Relative gap susceptibility as a function of electrode 
aspect ratio (horizontal axis: ratio of electrode thickness h to 
gap width g), substrate dielectric constant e s (line styles), and 
gap curvature (colors). Black: translation-invariant numeri- 
cal calculations of the setup of Fig. [3| with thick electrodes 
as in Fig. [fi] ao = —g/(2mv) [Eq. |14|]. For g <C w, a/ao 
does not depend on the strip width w. Colors: cylindrically 



symmetric numerical calculations of the setup of Fig. Ill with 

2 — 
Ri — and vanishing inner gap; qq = t§jt>. In : 



g 

32R ■ 



ponent similar to Eq. ([6| and a polarization potential 
equal to Eq. ^ , whose amplitude a is found from a line 
integral of the numerically determined electric potential 
across the gap. Comparing these numerical values of a 
to those known from analytic thin-electrode calculations 
(ao) yields the relative gap susceptibility a/ao- In com- 
plex setups this relative gap susceptibility cannot be ex- 
pected to depend only on local gap properties. But since 
Eq. ( 4b ) is dominated by a local quasi-singular contribu- 



tion, we expect the leading-order dependences of a/ao to 
be local, namely (i) the gap aspect ratio h/g, (ii) the sub- 
strate dielectric constant e s , and (iii) the gap curvature. 
Fig. [7] summarizes our numerical results showing these 
three dependences for several highly symmetric systems 
with uniform gaps. The main dependence is found to 
be on the gap aspect ratio. In the limit of large gap 
aspect ratios (h/g > 0.7) all of our results are in the 
vicinity of a/ao ~ 0-5, which leads us to propose the fol- 
lowing simplified model for gap simulations with signifi- 
cant aspect ratios: after numerically determining the gap 
polarizations ao by the thin-electrode inversion method 
of Section |II D| we simply multiply them by 0.5 before 
evaluating their electric potential. While this method 
does not take the exact gap curvature into account and 
incorrectly attributes equal weights to the influences of 
adjoining and far-off electrodes on the gap polarization 
through Eq. (4b), it nonetheless gives a good approxi- 
mation of the gap potential. Keeping in mind that the 
effects of gaps on three-dimensional potentials are gener- 
ally small (see Fig. [5]) , such an approximation is expected 
to be sufficient for the design phase of surface electrodes. 
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III. FINITE ELECTRODE PLANE 



notice that surface charge densities of the form 



The second approximation usually made for surface- 
electrode ion traps is that of an infinite grounded 
electrode plane surrounding the trapping region. But in 
realistic systems, the electrode plane is necessarily finite. 
In this section we attempt to take this fact into account 
for a particular electrode geometry in order to estimate 
its importance. 

Instead of assuming an infinite grounded plane we start 
from a planar grounded circular electrode of radius S 
surrounded by vacuum. All rf and dc electrodes will be 
placed on G. The reason why we choose to be circular 
is to simplify the following analysis in cylindrical coordi- 
nates; but substituting a rectangular for a circular one 
of similar size is not expected to change the conclusions of 
this section qualitatively and the tools developed in this 
section can thus be applied to a larger set of problems. 
For a given in-plane potential </>(r, $) on 0, determined 
by the shapes and potentials of the electrodes placed on 
0, the 3D potential $(r, z) must satisfy the following 
conditions: 

(i) The surface limit lim z _j. *&{ r , z ) — ( t ) ( r i $) must 
be valid for all r < S. However, for r > S we have 
no restrictions on the surface potential, as by as- 
sumption there is no electrode to fix the potential. 
The in-plane potential at r > S is to be determined. 



(ii) The surface charge density a(r, ■&) = 
— 2eq lim z _>o+ <9$(r, i9, z)/dz must vanish for 
r > S since by assumption there is no electrode to 
carry the charge density. But for r < S it is to be 
determined. 



As in Section II A there is a unique potential which sat- 
isfies both of these constraints. Here we derive this po- 
tential for an infinitesimally small "pixel" of nonzero in- 
plane potential on at radius p, given by 0(r, i?) = 
p~ 1 S(r — p)S(i}) with < p < S. The surface charge 
density is calculated from Eq. (4b). For r > p it can be 
multipole-expanded as [T7] 



a(r > p, $) 



E 

(m,n) 



7r(p 2 + r 2 — 2pr cost?) 3 / 2 

_> r ''' cos (mi?) 



2 3 -*™,o go ( n-rn-1 ) i ( n+m-1 ) % p n-3 



(22) 



where the sum is over m = 0, 1, 2, . . . and n = m + 3, m + 
5,m + 7, . . .; (<z)„ = T(a + n)/T(a) is the Pochhammcr 
symbol. This surface charge density must be brought 
to zero for r > S in order to satisfy condition (ii), but 
without modifying the surface potential for r < S. We 
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yield in-plane potentials through Eq. (4a) of 
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(24) 



where 2-F1 is Gauss' hypergeometric function and 3? refers 
to its real part. Adding such multipole charge densities 
to Eq. (22) thus allows us to null the charge density for 



r > S while leaving the associated surface potential on 
the electrode (r < S) unchanged. While we do not 
have a closed expression for the three-dimensional electric 
potentials due to the charge densities of Eq. (23), their 
series expansions for r 2 + z 2 < S 2 are 



$S,m,n(n z ) = " cos(toi?) 
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(25) 



Combining Eqs. ( 22 ) and ( 25 1 gives the three-dimensional 



potential of the in-plane potential "pixel" at radius p on 
the finite electrode 0, 
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cos(m^). (26) 



We recognize the first term in Eq. |26|) as the Green's 
function of Eq. ([2| in cylindrical coordinates, to which 
the 3D potential must converge in the limit S — > 00. 
The second term is the correction due to the finiteness 
of the electrode 0. The quadruple sum over Pochham- 



mer symbols in Eq. ( 26 ) can be reduced to a triple sum 



over hypergeometric functions by performing one of the 
summations analytically; however we have found that 
this does not make numerical evaluation faster in prac- 
tice. Arbitrary electrode patterns within can be con- 
structed from Eq. (JlJ by interpreting Eq. ( 26 1 as a mod- 
ified Green's function which takes the finite electrode 





FIG. 8: Perspective view of a planar ring ion (Paul) trap on 
a finite electrode. An rf electrode (green) is embedded in a 
finite grounded plane (gray) and surrounded by vacuum. 



A. ring ion trap on finite electrode 

As an example of a surface electrode pattern sur- 
rounded by a finite grounded electrode, we again study 
the optimized ring ion trap. Fig. [8] shows the schematic 
setup consisting of a ring-shaped rf electrode on a 
grounded disc of finite radius. Its three-dimensional elec- 
tric potential 

$s(r,z)= P dp [ MG s , p (r,#,z) (28) 

JRi JO 

can be simplified on the axis (r = 0) to 
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FIG. 9: Optimized planar ring trap with gapped electrodes, 
trapping at height z above the electrode plane. As a function 
of the reciprocal of the grounded electrode radius S, the opti- 
mal radii and the maximal achievable curvature k differ from 
the idealized result at S = oo. 



into account. The radius of convergence of Eq. (26 1 is 



r 2 + z 2 < S 2 , which is sufficient for practical simulations. 
For large S the Green's function in Cartesian coordinates 
is approximately 



G s (x,y,z) 



27r(x 2 + y 2 + z 2 ) 3 /2 3^25-3 



(27) 



adding to Eq. ^ a small constant out-of-plane electric 
field which can be easily incorporated into numerical cal- 
culations. 



As in Section II C Eq. ( 29 1 allows us to calculate the op- 
timal rf ring radii as a function of the size of 0: for given 
radius S and desired trapping height z we determine the 
radii R\ t 2 which maximize the dimcnsionless trap curva- 
ture k defined in Eq. ( 20 ) . Fig. [9] shows these quantities 
relative to their well-known values in the limit of an in- 
finite electrode plane S = 00. While the inner radius 
Ri is very insensitive to the finiteness of the substrate, 
the outer radius shows a rather strong dependence. This 
can be understood by the physical proximity of the outer 
radius to the substrate edge. As in Fig. [5j k is reduced 
from its ideal value (S — > 00) by only a few percent. 



IV. CONCLUSIONS 

After deriving general frameworks we have estimated 
the influence of finite gaps (Section II C ) and a finite sur- 



rounding electrode (Section III A) on the optimization 
of a ring-shaped Paul ion trap constructed from surface 
electrodes. These influences were found to be small in 
realistic situations: even extreme examples, such as a 
gap width of half the trapping height and a vanishing 
grounded electrode surrounding the rf trapping electrode, 
each give reductions in the dimensionless trap curvature 
of the order of only 5%. A similar small reduction was 



found in the far-field of a strip electrode (Section II B I. 
More realistic gap geometries including a finite gap as- 
pect ratio and a substrate dielectric (Fig. [6]) further re- 
duce the influence of gap potentials by a factor of two. 

We conclude that for many practical calculations of 
surface electrode designs it is a safe assumption to ne- 
glect the influence of both gaps and surroundings, since 
fabrication tolerances and stray fields will often intro- 
duce uncertainties in excess of these effects; practical 
implementations of surface ion traps necessitate control 



9 



electrodes for correcting these uncertainties, and can si- 
multaneously correct for design approximations such as 
those stemming from the effects quantified in this work. 

For those situations where higher accuracy is required, 
however, we have presented a method for calculating 
the effects of finite gaps and finite surrounding elec- 
trodes. This method provides fast estimates of the 
leading-order influences of gaps and finite electrodes on 
three-dimensional electric potentials generated by surface 
electrodes. The small size of these influences suggest that 
the two effects can likely be estimated separately, since 
the cross-polarizations of gaps and the "outside" plane 
are expected to be on the order of a few percent of the 
dominant polarization effects. After the design phase of a 



surface electrode setup using this approximate method, 
a full BEM calculation can then serve to validate the 
approximations and provide a fully accurate system de- 
scription. 
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